Constructing the phase diagram of finite neutral nuclear matter 
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The fragment yields from the multifragmentation of gold, 
lanthanum and krypton nuclei obtained by the EOS Collab- 
oration are examined in terms of Fisher's droplet formalism 
modified to account for Coulomb energy. The critical ex- 
ponents a and r and the surface energy coefficient cq are 
obtained. Estimates are made of the pressure-temperature 
and temperature-density coexistence curve of finite neutral 
nuclear matter as well as the location of the critical point. 



I. INTRODUCTION 

In past attempts to investigate the relationship be- 
tween nuclear multifragmentation and a liquid to vapor 
phase transition p]— [T^] various studies have sought to 
determine one or more critical exponents 0j3||] pT|^3| , 
other studies have examined caloric curves W, and still 
others have reported the observation of negative heat 
capacities |IJ. These studies suffer from the lack of 
knowledge of the system's location in pressure-density- 
temperature (p, p, T) space. For example, interpretations 
of caloric curves and negative heat capacities depend on 
assumptions of either constant pressure or constant den- 
sity In the case of determining critical exponents, 
it was assumed that the fragmenting system is at coex- 
istence and the dominant factor in fragment production 
was the surface energy. The analysis presented below 
makes no assumptions about the location of the system 
in (p, p, T) space and allows for other energetic consider- 
ations with regards to fragment production. 

In this paper the analysis technique recently used on 
multifragmentation data collected by the ISiS Collabora- 
tion is applied to the data sets for the multifragmen- 
tation of gold, lanthanum and krypton nuclei collected 
by the EOS Collaboration. All three EOS experimental 



data sets are shown to contain the signature of a liquid 
to vapor phase transition manifested by the scaling be- 
havior predicted by Fisher's droplet formalism, and the 
liquid-vapor coexistence line is determined over a large 
temperature interval extending up to and including the 
critical point. The critical exponents r and er, as well 
as the critical temperature T c , the surface energy coef- 
ficient Co and the compressibility factor Cf are directly 
extracted. From the behavior of the fragment yields the 
p-T and T-p coexistence curves are determined and the 
critical pressure p c and critical density p c are estimated. 



A. Overview 



The paper is organized as follows: section IB 



the EOS data sets; section II A reviews Fisher's droplet 



for malism and its connection to nuclear evaporation; sec- 
tion II B discusses the details of the data analysis; section 



II Oj reports the results of the data analysis; section |II D 
shows the physical implications of these results; and fi- 
nally, in section III a brief discussion of the results is 
made. In order to demonstrate the efficacy of the analy- 
sis performed on the EOS data sets, an appendix shows 
the results of this analysis performed on percolation clus- 
ter distributions. 



B. EOS data sets 

The EOS Collaboration has collected data for the re- 
verse kinematics reactions 1.0 AGeV Au+C, 1.0 AGeV 
La+C and 1.0 AGeV Kr+C There were 

-25,000, -22,000 and -36,000 fully reconstructed 
events recorded for the Au+C, La+C and Kr+C re- 
actions, respectively. The term "fully reconstructed" 
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means that the total measured charge in each event was 
within three units of the charge of the projectile. 

For every event, the charge and mass of the projectile 
remnant (Zq, Aq) were determined by subtracting the 
charge and mass of the particles knocked out of the pro- 
jectile from the charge and mass of the projectile p6| , p7t . 

The thermal component of the excitation energy per 
nucleon of the remnant E* was determined as follows. 
First, the total excitation energy per nucleon E^ otal was 
reconstructed based on an energy balance between the 
initial stage of the excited remnant and the final stage of 
the noninteracting fragments. The prescription [|l8| for 
calculating E? t , is then 



p* 

^total — 



£ {KEi + Q t ) + -nT 



/A 



(1) 



where n is the multiplicity of neutrons produced via frag- 
mentation, KEi is the kinetic energy of the i th fragment 
in the reference frame of the remnant and Qi is the re- 
moval energy and T is the temperature of a Fermi gas. 
For further details see reference Jl6[ . 

The thermal component of the excitation energy per 
nucleon of the remnant E* was then determined by sub- 
tracting the expansion energy Ex from E^ otal , where the 
quantity Ex is given by 



Ex — KE — £"th — ^Coulomb J /Aq 



(2) 



with KE the total kinetic energy; E t h the sum of the 
translational thermal contribution to the fragment spec- 
tra; and -Ecouiomb the Coulomb contribution. 
The translational energy is given by 



E t h = 2 T 'isotopc(w2 + n — 1) 



(3) 



where m,2 is the multiplicity of fragments and Ti sot0 po 
is the temperature calculated from the isotopic yields 
|l6| , [i"7f . This form follows that outlined in reference jl9) . 
The Coulomb contribution is given by 



E, 



Coulomb 



7 2 
R, 



E 



N,Z 2 



(V r /V )^R Z 



(4) 



where R r is the radius of the excited remnant, Nz is the 
number of fragments with charge Z, Rz is the radius (at 
normal density) of a fragment with charge Z. The vol- 
umes (and radii) were Vq the volume of the remnant at 
normal density and V r = Vp ro jcctiic(?Fer m i g a S /?lsotopo) 3/2 
the volume of the excited remnant that isentropically ex- 
pands from the normal volume of the projectile Vp ro jectiie 
p6[ ; R r is then determined from V r assuming a spherical 
volume. This form of -Ecouiomb follows reference 19 and 
takes into account the changing volume of the excited 
remnant as a function of excitation energy. Previous es- 
timates of -Ecouiomb did not account for the changing vol- 
ume of the fragmenting remnant [flq,[l7| . This difference 
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FIG. 1. (a) The distribution of events as a function of ex- 
citation energy, (b) The nucleon number of the fragment- 
ing system Ao normalized to the nucleon number of the pro- 
jectile ^projectile- (c) The charge of the fragmenting system 
Zo normalized to the charge of the projectile 

^projectile ■ (d) 

The number of intermediate mass fragments (4 < Z < Zo/A) 
maiF normalized to the charge of the fragmenting system, (c) 
The nucleon number of the largest fragment A ma x normalized 
to Aq. (f) The charge of the largest fragment Z max normal- 
ized to Zo. (g) The multiplicity of fragments (1 < Z < Z max ) 
resulting from the fragmentation of the system 7712 normalized 
to Z . 
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leads to a few AMeV difference in E* in the most violent 
collisions. 

For the analysis in this paper, the data for each system 
was binned in terms of E* in units of 0.5 AMeV; i.e. 
20 bins covered the excitation energy range OAMeV < 
E* < lOAMeV. Figure |l| shows some of the systematics 
of the EOS data binned in this manner. These results 
are consistent with other EOS publications p^Jl^ , [l7| , pO | . 

The systematics of the EOS data sets shown in Fig. l] 
demonstrates the similarity in behavior exhibited by the 
data sets when their differing sizes are taken into account 
by normalizing the quantity in question to the projectile 
charge ^projectile or the charge of the fragmenting system 
Zq- The exception is seen in Fig. |l|d where only below 
E* ~ 4 AMeV do all three systems behave similarly. 
Above E* ~ 4 AMeV, the size of the fragmenting systems 
dominate. This is reflected in the ordering of itiimf/Zo; 
from lowest to highest: krypton, lanthanum and gold. 



II. ANALYSIS 

As with several other analyses [0J|0-|ll|,||l| ||| , the ba- 
sis of the present effort lies in an examination of the frag- 
ment yield distribution in the context of Fisher's droplet 
formalism [p5|-p9|. Thus, a brief review of Fisher's for- 
malism is given in the following section, together with a 
justification for its applicability to nuclear decay rates. 

A. Fisher's droplet formalism 

Fisher's droplet formalism and its forerunners [ p0|j3l] 
are based on an equilibrium description of physical clus- 
ters or droplets that condense in a low density vapor. 
While Fisher's formalism has long been applied to nuclear 
multifragmentation yields []l1,^|J^- [Tl| , |2l| -p4|] , the question 
arises as to the validity of a picture of clusters in equilib- 
rium within a low density vapor to experiments in which 
excited nuclei undergo multifragmentation in vacuum. 
Specifically: "In which sense is there an equilibrium be- 
tween liquid and vapor in the free (vacuum) decay of 
a (multifragmenting) hot intermediate (nucleus)?" Or 
more to the point: "Where is the vapor?" 

If one assumes, as in a compound nucleus reaction, 
that the initial collision entity relaxes quickly to a hot 
thermalized blob, which proceeds slowly to emit particles 
stochastically, the answers to this question is simple. The 
hot blob is the liquid which is evaporating in free space 
according to standard evaporation theories. To establish 
coexistence, the vapor need not be present. All that is 
necessary is to appreciate that: 

1. in first order phase transitions the interaction be- 
tween the two phases is unnecessary; 



2. the rate of evaporation defines uniquely the vapor 
phase, even when the vapor phase is absent. 

In fact the concentration of species A is completely de- 
fined by 



Ra = n A (T)v A {T) 



(5) 



where Ra(T) is the emission flux of A, ua(T) is the con- 
centration of species A and WX{T) is the average velocity 
of A which is of order of yjT/A. In other words, the 
outward flux, at equilibrium, is the same as the inward 
flux. 

Thus a direct connection is made between the statis- 
tical decay rate and Fisher's equilibrium description of 
cluster formation. Two consequences follow: 

1. At equilibrium, the evaporated particle is replaced 
by the back flux from the vapor. However, since the 
back flux is absent in the case of nuclear multifrag- 
mentation, this analysis is limited to particles with 
low emission probability (first chance) and must 
avoid particles which are emitted with high mul- 
tiplicity. This is approximately achieved by elim- 
inating fragments with Z < 4 from the ensuing 
analysis. 

2. In the same spirit as above, the pertinent tempera- 
ture is that of the blob as it evaporates low proba- 
bility particles. Thus, rather than worrying about 
the role of high multiplicity particles and their asso- 
ciated cooling on the energy-temperature relation- 
ship, the Fermi gas relation ship E — aT 2 can be 
assumed with good confidence. 

With this picture in mind, we return to Fisher's formal- 
ism. 

The basic idea is that in a non-ideal vapor of particles 
interacting with repulsive cores and short range attrac- 
tive forces, can be approximated at low densities and 
temperature by an ideal gas consisting of noninteracting 
monomers, dimers, trimers at equilibrium. The (free) 
energy of sufficiently large clusters can be estimated in 
terms of their volume and surface energy. These clus- 
ters are in equilibrium with each other and the relative 
abundances of differently sized clusters changes with tem- 
perature and pressure |2(| . 

The relative abundances of clusters with A con- 
stituents is given by: 



n A = q A T exp 



AfiA c eA° 



T 



T 



(6) 



Here qo is a constant of proportionality which is fixed 
by the critical density [ p2p^ ]. The power law A~ T arises 
from a combinatorial factor that depends on the fact that 
the surface of the cluster must be closed | p4p5| . The 
distance from coexistence is: 
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(7) 



where /i/ is the chemical potential the liquid at coexis- 
tence and fi is the chemical potential of the system. For 
A/i > (a super-heated vapor) and A/x = (liquid- vapor 
coexistence) the above sum always converges. While for 
A/i < (a super-saturated vapor), the sum diverges. 
The "classical" part of the surface energy is parameter- 
ized by cqsA' 7 , where cq is the zero temperature surface 
energy coefficient, s — (T c — T)/T c and A a relates the 
number of constituents of a cluster to the most probable 
surface area. Fisher's critical exponents a and r depend 
on the Euclidean dimensionality and universality class of 
the system. 

The total pressure of the entire cluster distribution is 
given by summing all of the partial pressures Tua- 



P 



= zZ TnA 



and the density is 



p = y^AriA- 



(8) 



(9) 



Thus the pressure and density of the system can be in- 
ferred from the knowledge of the cluster distributions. 

At the critical point the system is at coexistence (A/i = 
0) and the classical part of the surface energy cost van- 
ishes (e = 0). Thus both exponential factors are unity 
leaving only the temperature independent power law: 



n A = qoA T . 



(10) 



Away from the critical point, but along the coexistence 
curve so that A/i = 0, the cluster distribution is given 
by: 



n A = q A T exp 



c eA° 
T 



(11) 



Equation ( |ll| ) can be rewritten as: 

n A = QoA exp — — exp 



T 



i?exp 



(12) 



Thus the cluster distribution along the coexistence curve 
is given by a Boltzmann factor with 



fc A° 
K = qoA exp ' 



and 



B = c A". 



(13) 



(14) 



This Boltzmann factor manifests itself in Arrhenius plots 
for the fragment yields where a linear relation between 



ln^^) and and 1/T is observed. This behavior has long 
been observed in many nuclear fragment yield distribu- 
tions and has recently been observed in the clus- 
ter distributions of percolation (with the bond breaking 
probability playing the role of temperature) ]3(| and the 
Ising model Q . 

As discussed above, Fisher's formalism relates directly 
to a reaction rate picture. In this picture, the heavy frag- 
ments (e.g. Z > 4) are the product of first chance emis- 
sion from the excited remnant. The first chance emission 
from a compound nucleus can be written as: 



n^Tjocrae^. 



(15) 



Thus the fragment yields, parameterized via Fisher, can 
be related to the decay rates (widths T). Furthermore, 
r, which controls the first chance emission yields, is the 
same decay width which controls the mean emission times 
t since: 



thus 



Tt 



t OC — (X — - 

r n A (T) 



,(B/T) 



(16) 



(17) 



The mean time for fragment emission reported by the 
ISiS Collaboration [ ^l|j43| is well described as a Boltz- 
mann factor. It was also noted that the Boltzmann fac- 
tor describing the emission times is the same as that de- 
scribing the fragment yields (44| . This indicates that the 
thermal reaction rate picture is valid for multifragmenta- 
tion; fragments can be viewed as being the result of the 
evaporation of an excited nucleus. 



B. Fitting the data 

Preliminary fits of the gold, lanthanum and krypton 
data with Eq. (^]) led to puzzlingly large results for A/i 
((A/i) w 3 AMeV) which could be interpreted as a sub- 
stantial degree of super-saturation. A much more plau- 
sible alternative explanation is the lack of account of the 
Coulomb effects in Fisher's formalism. Equations ( fl2] ) 
and (|l4j) support the presence of a barrier controlling 
the flux from liquid to vapor and vice versa. This bar- 
rier should depend not only on the surface energy of the 
fragment but should reflect the entire energy necessary 
to remove a fragment from the liquid and place it into 
the vapor. At the least, the energy necessary to relo- 
cate a charge Z from the bulk to "near" the surface of 
the "residual" nucleus should be evaluated. This energy 
is negative and counteracts the effects of the surface en- 
ergy. Furthermore, it is to leading order, linear in Z and 
thus in A. Thus the large values of A/i. An attempt to 
include the Coulomb energy explicitly is then 
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n A = q A T exp 
with 

-^Coul 



AA/j + Ecoui c eA a 



T 



(Zo -Z)Z 



T 



(18) 



r ((A - A) 1/3 + AV3) 



(1-e— ) (19) 



with ro = 1.2 fm. In Eq. (19) there is a recognizable 
Coulomb interaction energy of two touching spheres mod- 
ified by a factor (1 — exp(— xe)). The parameter x (left 
as a fit parameter) takes into account the numerical co- 
efficients of the linear term in Z plus polarization effects, 
and e takes care of the need for the vanishing difference 
between the liquid and vapor near the critical point. Note 
that the Coulomb energy discussed in Eq. (||) is different 
from the Coulomb energy discussed in Eq. (|l9|). Equa- 
tion (Q) describes the total Coulomb energy present in 
the fragmentation process, while Eq. (|l^) describes the 
cost in moving a fragment from the nuclear liquid to the 
nuclear vapor. 

The mass of a fragment A prior to secondary decay was 
estimated by multiplying the measured fragment charge 
Z by 2 and then by a factor of (1 + y(E*/Bf)) where Bf 
is the binding energy of the fragment and y is a fit param- 
eter to allow for an increase or decrease in the amount of 
secondary decay. 

The temperature was determined by assuming a de- 
generate Fermi gas, 



T = VaE*. 
The parameter a was taken to be [M5| 



a 



1 



E* 
So 



(20) 



(21) 



in order to accommodate the empirically observed change 
in a with excitation energy [ fl6|| . Here Bq is the binding 
energy of the fragmenting system. Using the Fermi gas 
approximation to relate E* and T gives a reasonable es- 
timate of the temperature of the excited remnant at the 
time of first chance emission. It has been observed that 
even the isotope ratio thermometer follows the Fermi gas 
approximation quite well so long as the average number 
of IMFs is less than one jl7) . 

To obtain the concentration of fragments of a given 
mass, the total number of fragments Na of a given size 
A was normalized to the size of the fragmenting system 
A so that ha — Na/Ao. 

The location of the critical point, in terms of exci- 
tation energy, was determined from an examination of 
measured fluctuations. In general, as the critical point 
of a system is approached from the two phase region, 
the difference between phases diminishes and the system 
fluctuates from one phase to the other. At the critical 
point the fluctuations are maximal. However, while the 



maximum in the fluctuations occurs at the critical point, 
the presence of a peak in the fluctuations is a necessary, 
but not sufficient, condition for a existence of a phase 
transition 0. 

The fluctuations measured in the EOS data are: (1) 
in the charge of the largest fragment normalized to the 
charge of the fragmenting system, and (2) related to the 
average mass number of a fragment as measured by the 
quantity 72 [E8| , where 



72 



/RMS(A) \ 2 + i _ M 2 M 



(22) 



V <^> J ' * Ml 
with Mfc as the fc th moments of the fragment distributions 



AI k = ^n A A k . 

A=l 



(23) 



These fluctuations are shown in Fig. ||. 

The peak in the fluctuations was found by smooth- 
ing the data (solid lines in Fig. |^), taking the numerical 
derivative of the smoothed data and finding the value 
of E* where the the derivative passed through zero, see 
Fig. ||. Finally, the value of the excitation energy at the 
critical point E* was determined by averaging the results 
from both measures of the fluctuations. Table | lists the 
results. For this analysis the values determined for the 
excitation energy at the critical point for the Au+C reac- 
tion are in proximity of other values observed in previous 
EOS analyses (E* » 4.75 AMeV) |§||||l7). Differences 
in the values of E* arise from the different method of 
constructing the thermal portion of the excitation energy 
described above fl2| , po|| . 

Estimates of the critical temperature T c are made by 
using the values of E* in Eq. ( pp| ) and lead to values, 
shown in Table Q, that are comparable to theoretical es- 
timates for small nuclear systems p9|-p2]|. As an aside, as 
shown in Table | the value of T c increases with decreasing 
projectile (and thus remnant) mass. This is opposite of 
the trend assumed in a prior analysis of the EOS gold 
multifragmentation data where the Coulomb energy was 
neglected || but in agreement with the trend reported 
in other analysis of the EOS data sets |l2| , p0| ]. 

In Fig. |^ the value of 72 for the Kr system attains a 
peak value of only ~1.8. It has been suggested that the 
magnitude of the peak in 72 could distinguish between 
the presence of a power law with r > 2 (72 > 2) and 
an exponential distribution (72 < 2) in the cluster yields 
|^,^8|. However, it was seen that this is not the case || 
and it will be seen in the Appendix that small percolation 
lattices have values of 72 with peak magnitudes of less 
than two yet still exhibit a continuous phase transition 
with an exponent of t <~ 2.2 in the power law describing 
the cluster yields at the critical point. Thus, the height of 
the peak in 72 cannot be used to rule out the presence of a 
critical point and the associated power law in the cluster 
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E* (AMeV) E* (AMeV) 

FIG. 2. Left: The RMS fluctuations in the charge of the 
largest fragment normalized to the charge of the fragmenting 
system plotted as a function of excitation energy. Right: The 
quantity 72 plotted as a function of the excitation energy. 
Open symbols show the data points, solid curves show the 
results of smoothing the data. 



distribution or provide information about the value of 
the power law exponent. 

Data from each system for 0.25 AMeV < E* < E* 
(which corresponds to a range of < e <~ 0.8) and 
4 < Z < Zq/4 were simultaneously fit to Eq. ([i"8|), which, 
as mentioned previously, helps insure that the fragments 
examined in this analysis are produced via first chance 
emission. There were nearly 200 points from the EOS 
data sets used in the fitting procedure. The fit param- 
eters r, a and cq were kept the same for all three data 
sets while Ap:, x and y were allowed to vary between 
the systems to minimize chi-squared; this gives 12 free 

TABLE I. Critical points of excited nuclei 



System 


E* (AMeV) 


T c (MeV) 


pc (po) 


p c (MeV/fm 3 ) 


Au + C 


4.6 ±0.2 


7.6 ±0.2 


0.39 ±0.01 


0.11 ±0.04 


La + C 


4.9 ±0.2 


7.8 ±0.2 


0.39 ±0.01 


0.12 ±0.04 


Kr + C 


5.1 ±0.2 


8.1 ± 0.2 


0.39 ±0.01 


0.12 ±0.04 




2.5 5 7.5 10 2.5 5 7.5 10 

<£*> (AMeV) <£*) (AMeV) 

FIG. 3. Numerical derivatives of the smoothed data from 
Fig. ^ of: (left) the RMS fluctuations in the charge of the 
largest fragment normalized to the charge of the fragmenting 
system plotted as a function of excitation energy, and (right) 
the quantity 72 plotted as a function of the excitation energy. 
Solid star symbols shows the approximate location in excita- 
tion energy where the derivative is zero, thus indicating the 
critical point. 



parameters used to fit nearly 200 data points. Previous 
analyses of the EOS data |||36| assumed that A/i = 
and that the effects of the Coulomb energy were small. 
The analysis presented here makes no such assumptions. 

Fixing t at 2.2 did not significantly change the results 
of this analysis. Using a common x value for all three 
data sets also returned results similar to those quoted be- 
low. Using a common y value for all three data sets also 
returned results similar to those quoted below. These 
different methods suggest a systematic error of ~ 15% of 
the value in question. All errors quoted below are those 
returned by the fitting procedure, propagated where nec- 
essary. Finally, the same data collapse observed below 
would be seen if the parameters were fixed to: r = 2.21, 
a = 0.64 (their d = 3 Ising values), c = 16.8 MeV 
(the text book value of the nuclear liquid-drop surface 
energy coefficient), A/i Au LaKr = (they must be close 
to zero since fragments are observed) and y — 0.5 (in 
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TABLE II. Uncommon fit parameters 



System 


Afi (AMeV) 


X 


y 


Au + C 


0.38 ±0.02 


1.1 ±0.2 


0.5 ±0.1 


La + C 


0.47 ±0.03 


1.2 ±0.1 


0.3 ±0.2 


Kr + C 


0.58 ±0.08 


4.0 ± 1.0 


0.8 ±0.2 



keeping with previous assumptions that the fragments 
prior to secondary decay have the same mass to charge 
ratio of the excited remnant P,^2[) and letting only x, 
the Coulomb parameter vary to minimize chi-squared. 



C. Results 

Figure |i] shows the fragment mass yield distribution 
scaled by the power law pre-factor, the chemical potential 
and Coulomb terms: riA/qoA~ r exp((A^v4 + Ecou\)/T) 
plotted against the inverse temperature scaled by Fisher's 
parameterization of the surface energy: A a e/T. Now, 
the scaled data for all three systems collapse onto a sin- 
gle line over several orders of magnitude as predicted by 
Fisher's droplet formalism p5|] . This collapse provides 
direct evidence for a liquid to vapor phase transition in 
excited nuclei. Furthermore, the fact that the data from 
each system show a common scaling illustrates the com- 
mon nature of the underlying phenomenon. 

1. Parameters 

The values of r = 2.2 ± 0.1, a = 0.71 ± 0.02 and 
co = 14.0 ± 1.0 MeV determined in this analysis are in 
agreement with those determined for the ISiS gold mul- 
tifragmentation data sets fl"c|| and are in agreement with 
values previously determined for the EOS Au±C data set 
[^2|,^6|. The value of the surface energy coefficient cq is 
close to the value of the surface energy coefficient of the 
liquid-drop model which is ~ 16.8 MeV. 

A previous analysis of the EOS gold multifragmenta- 
tion showed the surface energy coefficient to be cq = 6.8± 
0.5 MeV @. The difference between the c = 6.8 ± 0.5 
MeV from that work and the c = 14.0 ±1.0 MeV pre- 
sented here arises from the differing analyses. In the 
previous analysis it was assumed that A/x = 0, that the 
Coulomb energy was negligible and that the level density 
parameter was constant at a = 13. These assumptions 
allowed some degree of scaling, yielded sensible values for 
the critical exponents but resulted in a surface energy co- 
efficient that was a factor of two of lower than that of the 
present analysis. 

In addition to a surface energy coefficient that is in bet- 
ter agreement with the standard liquid-drop model, the 
greater collapse of the data in the present work demon- 
strates the improvements of the present analysis over the 



previous one. The improvements in analysis are related 
to allowing a non-zero A/x, taking into account the cost 
in Coulomb energy to move a fragment from the liquid 
to the vapor and accounting for the change in the level 
density parameter over the excitation energy range. The 
treatment of secondary decay in both analyses is differ- 
ent: previously it was assumed that the fragments, prior 
to any secondary decay, had the same mass to charge ra- 
tio as the fragmenting remnant. In the present analysis 
the amount of secondary decay is left as a free parameter. 

The values of A/x reported in Table || can be consid- 
ered "small" in light of Eq. . The chemical potential 
of the liquid can be found by 



Hi = E + TS 



(24) 



with Eq as the bulk energy per particle and So as the 
bulk entropy per particle |25[ . Treating the system as a 
Fermi gas so that S/A = aT yields 



/x = E + E* 



(25) 



Thus the important energy scale for A/i is Eo+E*, for nu- 
clear matter Eq ~ 15.5MeV. The values of A/i returned 
by this analysis are < 6% of Eq ± E* indicating the sys- 
tem is close to coexistence. The values of A/x should also 
be compared to the values returned when the EOS frag- 
ment yields were fit to (Eq. ||): (A/i) w 3.0 AMeV for all 
EOS reactions. The reduction in the magnitude of the 
A/x values is about a factor of six and is due to the mod- 
ification of Eq. (Q) to account for the Coulomb energy, 
i.e. Eq. (|l8|). The remaining small positive A/x values of 
the systems may indicate that those systems are slightly 
super-saturated, or more probably they may reflect some 
other energy costs not taken into account (e.g. the sym- 
metry energy or pairing), or they may reflect that the 
approximation for the cost in Coulomb energy to form a 
fragment given in Eq. (|l^) is not completely adequate (for 
instance Eq. ( |l9|) assumes a spherical geometry which 
may or may not be the case), or they may merely reflect 
noise in the data. 

The values of x for each system may indicate more (Au 
and La) or less (Kr) Coulomb energy present in the sys- 
tem. They may also reflect the symmetry of the collision 
which may affect the geometry of the remnant, e.g. a 
very asymmetric collisions like Au±C may leave a nearly 
spherical remnant, while a more symmetric collision like 
Kr±C may result in a less spherical fragmenting system. 

The values of y returned indicate that the fragments 
have the same mass to charge ratio as the excited rem- 
nant. 

The difference in values of A/x, x and y determined in 
the analysis of the three EOS data sets and those deter- 
mined in the analysis of the ISiS 8.0 GeV/c it on gold 
multifragmentation set |10| is left an open question. The 
small differences in E* and T r are due to the differences 
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in reconstructed excitation energy scales |pq |. This dif- 
ference carries over to all energy related quantities, e.g. 
Co- 

Finally, in light of the above parameter results, it is 
clear that the same data collapse would be observed if the 
parameters were fixed to some nominal values, discussed 
above, with only x, the Coulomb parameter varying to 
minimize chi-squared. Thus only three free parameters 
are truly needed to fit the ~ 200 data points of the EOS 
data sets. 



D. The coexistence curve of finite neutral nuclear 
matter 

1. The pressure-temperature coexistence line 

Before determining the pressure-temperature coexis- 
tence line, the meaning of a pressure associated with an 
excited nuclear remnant must be addressed. As discussed 
above, in the actual experiment, this pressure is virtual; 
it is the pressure the vapor would have in order to pro- 
vide the back flow needed to keep the source at equilib- 
rium. However, since the yields from Fisher's formalism 
are proportional to both the pressure, Eq. (||), and the 
evaporation rate, Eq. (|l7|), it is clear that by fitting the 
yields as has been done above, one can infer an associated 
(virtual) vapor pressure. 

The p-T coexistence curve can be determined from this 
analysis. As seen in section II A, Fisher's theory assumes 



a, 



that the non-ideal fluid can be approximated by an ideal 
gas of clusters. Accordingly, the quantity ua is propor- 
tional to the partial pressure of a fragment of mass A 
and the total pressure due to all of the fragments is the 
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FIG. 5. The reduced pressure-temperature phase diagram: 
the points show calculations performed at the excitation en- 
ergies below the critical point and the lines show fits to the 
Clausius-Clapeyron equation. 



sum of their partial pressures (see Eq. 
pressure is then given by: 

p _ TJ2n A (T) 
Pc T c Y,n A (T c )' 



)). The reduced 



(26) 



The coexistence curve for finite neutral nuclear matter 
is obtained by substituting the Ua{T, A/i = 0,Ec O ui = 
0) from Eq. (|l8|) in the numerator of Eq. ( |26| ) and 
ua(T c , A/i = 0, Ecoui = 0) in the denominator. This al- 
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lows one to transform the information in Fig. || into the 
familiar phase diagram in Fig. ^. The data points shown 
give the values of p/p c and T c /T calculated via Eq. ( p6| ) 
for the bins in E* up to and including the critical point. 

Figure [5] gives an estimate of the coexistence line of 
finite nuclear matter and from this it is possible to make 
an estimate of the bulk binding energy of nuclear matter. 
One begins by assuming the system behaves as an ideal 
gas and uses the Clausius-Clapeyron equation 



dp_ _ AH 
dT ~ TAV 



(27) 



where AH is the molar enthalpy of evaporation and AV 
is the molar volume difference between the two phases. 
Then solving for the vapor pressure with 



AV = V, 



vapor 



-Vi 



iquid 



vapor 



T 
P 



gives 



-AH 

p = p exp(^^) 



which would lead to the ratio of 



P_ 

Pc 



exp 



AH 

~t7 



1 - 



T 



(28) 



(29) 



(30) 



if AH were assumed to be temperature independent. 
However, as T — > T c the gas is not ideal and AH ^ 
constant, but it has long been known that for several nor- 
mal fluids these deviations compensate so that ln(p/p c ) 
is approximately linear in T/T c fl54| . 

A fit of Eq. ( |30| ) to the coexistence curves for the sys- 
tems is shown in Fig.|5| yields the ratio of AH/T C . Using 
the corresponding values of T c gives the molar enthalpies 
of evaporation of the liquid AH shown in Table III . From 
these AH values AE is constructed via AE — AH — pV 
with pV = T (with the ideal gas approximation) using 
the average temperature from the range in Fig. I listed m 
Table [II. AE refers to the cost in energy to evaporate a 
single fragment. To determine the energy cost on a per 
nucleon basis AE is divided by the most probable size of 
a fragment over the temperature range in Fig. |^. Since 
the gas described by Fisher's formalism is an ideal gas of 
clusters, the most probable cluster size is greater in size 
than a monomer. The most probable size of a fragment 
in the region of the p-T coexistence line obtained from 
Eq. (^) and the experimentally determined parameters 
is 1.05 ± 0.05. Thus the AE/A becomes wl4 AMeV, 
close to the nuclear bulk energy coefficient of 15.5 MeV. 



TABLE III. Thermodynamic properties of excited nuclei 
System AH (MeV) (T) (MeV) AE/A (AMeV) Cf 



Au + C 19.4 ± 0.7 
La + C 19.6 ± 0.7 
Kr + C 19.5 ±0.7 



4.6 ± 0.6 14 ± 1 
4.9 ±0.6 14 ±1 
4.9 ±0.6 14 ±1 



0.28 ± 0.09 
0.28 ±0.09 
0.28 ± 0.09 



6 




FIG. 6. The points are calculations performed at the exci- 
tation energies below the critical point and the lines are a fit 
to and reflection of Guggenheim's equation. 



p_ = J2An A (T) 



(31) 



With Ap and Ecoul set to in the numerator of Eq. JlS| ) 
and Ap and Ecoul set to with T set to T c in the de- 
nominator, Eq. (|3l| ) gives the low density (vapor) branch 
of the coexistence curve of finite nuclear matter, shown 
in Fig. |. 

Following Guggenheim's work with simple fluids, it is 
possible to determine the high density (liquid) branch 
as well: empirically, the p/ p c -T/T c coexistence curves of 
several fluids can be fit with the function fp5|: 



Pl,v 
Pc 



= 1 + 61(1 



T T 



(32) 



where the parameter 62 is positive (negative) for the liq- 
uid pi (vapor p v ) branch. Using Fisher's formalism, /3 
can be determined from r and a [E5fl : 



(3 = 



(33) 



2. The temperature- density coexistence curve 



As seen in section [I A the system's density can be 
found from Eq. (j|). The reduced density is given by: 



For this work j3 = 0.3 ± 0.1. Using this value of j3 and 
fitting the coexistence curve from the EOS data sets with 
Eq. ([32|) one obtains estimates of the p v branch of the 
coexistence curve and changing the sign of 62 gives the 
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pi branch, thus yielding the full T-p coexistence curve of 
finite nuclear matter. 

From Fig. |6| it is possible to make an estimate of the 
density at the critical point p c . Assuming that normal 
nuclei exist at the T = point of the pi branch of the 
coexistence curve, then using the parameterization of the 
coexistence curve in Eq. ( |32| ) gives p c ~ po/3. See Table | 
for precise values. 

3. The pressure-density coexistence curve 




FIG. 7. The points are calculations performed at the exci- 
tation energies below the critical point and the lines are the 
results of the fits from the previous sections. 



For the sake of completeness the p/p c ~p/Pc projection 
of the coexistence curve is determined by combining the 
results of the previous two sections. This is shown in 
Fig. 0. It is clear from Fig. that the fitted curves do 
not reach p/p c = 1 at p/p c = 1 while the data points do. 
This is a reflection of the validity of the assumptions that 
went into deriving Eq. (|30"|). 

4- The compressibility factor 



Finally, a measure of the pressure at the critical point 
p c can be made by using T c and p c from above in combi- 
nation with . The results are shown in Table [jj. This 
last calculation gives a complete experimental measure 
of the location of the critical point of finite neutral nu- 
clear matter (p c , T c , p c ) and is in agreement with the ISiS 
results and in rough agreement with theoretical calcula- 
tions nm. 




0.5 0.6 0.7 0.8 0.5 0.6 0.7 0.8 

q q 
FIG. 8. Left: The RMS fluctuations in the size of the 
largest cluster normalized to the size of the lattice system 
plotted as a function of bond breaking probability. Right: 
The quantity 72 plotted as a function of bond breaking prob- 
ability. Open symbols show the estimate of the excitation at 
the critical point based on the maximum of the fluctuations, 
the solid line shows the results of smoothing the data. 



The critical compressibility factor = p c /T c p c can 
also be determined in a straightforward manner from p8[ : 



C* = 



EAn A (T c 



(34) 



Table III shows the results for the EOS data sets which 
are in agreement with the values for several fluids |2|] 
and that of the ISiS data ifLOl. 



III. CONCLUSION 

Through a direct examination of the most accessible 
features of nuclear multifragmentation, namely the frag- 
ment distributions themselves, and the use of Fisher's 
droplet formalism, modified to account for the Coulomb 
energy cluster formation, a measurement of the coex- 
istence curve of finite neutral nuclear matter has been 



10 




0.5 0.6 0.7 0.8 0.5 0.6 0.7 0.8 



<q) (q) 
FIG. 9. Left: The numerical derivative of RMS fluctua- 
tions in the size of the largest cluster normalized to the size 
of the lattice system plotted as a function of bond breaking 
probability. Right: The numerical derivative of the quantity 
72 plotted as a function of bond breaking probability. Open 
symbols show the estimate of the excitation at the critical 
point based on the maximum of the fluctuations, the solid 
stars show where the derivatives are zero. 

made for three different multifragmenting systems and 
estimates of the critical point for finite nuclear matter 
have been made. The precise values of quantities like the 
critical exponents and critical temperature and precise 
locations of coexistence curves depend on the assump- 
tions made for the cost in Coulomb energy for fragment 
formation and the assumptions made to account for the 
secondary decay of the fragments. While the exact forms 
are unknown, the estimates made in this paper have solid 
physical origins, and yield values of the surface energy co- 
efficient and the bulk binding energy of nuclear matter 
which are consistent with established values. Both the 
p-T coexistence lines and the T ~ p coexistence curves for 
all three EOS systems are consistent. These are strong 
indications that this analysis determines the coexistence 
curve and can be used to construct the phase diagram of 
finite neutral nuclear matter based on experimental data. 



TABLE IV. Critical points of finite percolation lattices 



L 


<7c 


pc 




9 


0.705 ± 0.004 


0.210 ±0.001 


0.041 ±0.001 


6 


0.685 ± 0.004 


0.216 ±0.001 


0.041 ±0.001 


4 


0.655 ± 0.004 


0.243 ± 0.002 


0.044 ± 0.001 



IV. APPENDIX 

To demonstrate the efficacy of the above analysis, it 
is applied to the cluster distributions from three dimen- 
sional simple cubic lattices of side L = 4, 6 and 9. It will 
be seen that if the above procedures are followed, well 
known quantities are recovered. 

Cluster distributions for over 100, 000 lattice realiza- 
tions were generated by breaking bonds between sites 
p6[ . A value of the lattice's bond breaking probability q 
was chosen from a uniform distribution on (0,1). Next, 
a bond probability, qt, was randomly chosen from a uni- 
form distribution on (0,1) for the i th bond. If qi was less 
than q, then the i th bond was broken and two sites were 
separated. This process was performed for each bond in 
the lattice. At low values of q, few bonds were broken 
resulting in a cluster distributions that are analogous to 
the liquid-vapor coexistence of a fluid. In an infinite lat- 
tice the distinguishablity of the "liquid" phase and the 
"vapor" phase vanishes at a unique value of the lattice 
probability, q c , when the probability of forming a perco- 
lating cluster changes from zero to unity J57],Q . For the 
ensuing analysis, the number of clusters of size A per lat- 
tice site ha was calculated by histogramming the lattice 
realizations into 100 bins on q from to 1. 

First the value of the probability at the critical point 
q c is determined by locating the maximum in the fluc- 
tuations of: (1) the size of the largest cluster, and (2) 
72. Figures || and |^ show these measures of the fluctu- 
ations. The location of the maximum is determined as 
in the EOS data, the data is smoothed and then the nu- 
merical derivative is taken. The location of the peak in 
the largest cluster is averaged with the location of the 
peak in 72 and the results are recorded in Table As 
expected the value of q c changes with the lattice size. 

Note that in Fig. || the value of 72 for the L — 4 lattice 
attains a peak value of only ~1.9; this is a finite size effect 
and due to the small size of the lattice. Since 72 is related 
to the fluctuations in the average size of a cluster, it is 
clear that as the size of the lattice decreases, the upper 
limit in the size of a cluster decreases, thus imposing a 
limit on the size of 72. 

Next the cluster yields from the three different lat- 
tices are fit simultaneously to Eq. (|J), with q c {L) keep- 
ing the fit parameters tr and r consistent between lattices 
and letting A/i and cq vary between lattices. Data from 
0.4 < q < 1.05q c and 5 < A < 3L were included in the 
fitting procedure. This gives seven fit parameters with 
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FIG. 10. The scaled yield distribution versus the scaled bond breaking probability for the L = 9, 6 and 4 lattices. The solid 
lines has a slope of co(L). 



TABLE V. Percolation fit parameters 



L 


A/i 


CO 


9 


-0.008 ± 0.004 


2.62 ±0.04 


6 


0.001 ± 0.001 


2.42 ± 0.04 


4 


0.007 ± 0.001 


1.91 ±0.04 



TABLE VI. "Thermodynamic" properties of finite perco- 
lation lattices 

L AH AE/A Cf bonds/site 

9 3.62 ±0.03 2.7 ± 0.1 0.275 ± 0.003 2J7 

6 3.35 ±0.03 2.2 ±0.1 0.275 ± 0.003 2.50 

4 2.75 ±0.03 1.8 ±0.1 0.275 ± 0.003 2.25 



1083 points to fit. The results are shown in Fig. [LO] and 
recorded in Table [v|. 

The formula in Eq. (J^j used in this analysis is only one 
example of a more general form of the scaling assumption 



n A =A- T f(X) 



(35) 



with X = A a e v /T and where f(X) is some general scal- 
ing function. This scaling function should be valid on 
both sides of the critical point. For small X (T w T c 
and small A) and e > 0, f(X) will vary as exp(— X) with 
cr = l/(/3S) = 1/(7 + (3) = 0.64 for d = 3 Ising systems 
or 0.45 for d — 3 percolation systems and (p = 1. For 
large X (T far from T c or large A) and e > 0, f(X) will 
vary as exp(— X) with a = 2/3 for all three dimensional 
systems and with tp — 1v\ where v — 0.63 for d = 3 Ising 
systems and v = 0.88 for d — 3 percolation lattices. 

The fitting procedure using Eq. (|^) returned a value of 
(7 = 0.44 ±0.01 and r = 2.192 ±0.003 in good agreement 
with other measurements, a — 0.45 and r = 2.18 |5^| , It 
is clear from these results that the data examined here is 
in the small X , s > region where the approximation of 
f(X) given in Eq. (||) is valid. As with the EOS data, the 
errors quoted here are from the fitting procedure. Sys- 
tematic errors that arise from the use of Fisher's scaling 



form and from the fitting regions in A and q are on the 
order of ~ 10%. 

The value of Co for the L = 6 lattice is in good agree- 
ment with previous measures fl36| . The interpretation of 
the change in cq with lattice size will be discussed below. 

The values of A/i for all lattices are close to zero, 
in agreement with the fact that percolation calculations 
such as these are at coexistence. 

It is now a simple matter to follow the analysis de- 
scribed above using Fisher's parameterization of the clus- 
ter distribution to determine the "phase diagrams" for 
these percolation lattices. The interpretation of these 
"phase diagrams" is not as simple. 

First the "reduced pressure" as a function of the in- 
verse of the "reduced probability" q/q c is determined via 
Eq. (prjj), and as usual for percolation studies q replaces 
T and q c replaces T c . The results are shown in Fig. |ll| 
where the points are fit with Eq (|3(]) . This leads to an es- 
timate of the "enthalpy of evaporation of a cluster" given 
in Table VI. The values of AH are on the order of the 



values of Co and increase with increasing L. 

To determine the "energy of vaporization" of a cluster 
AE the ideal gas approximation pV — q is followed so 
that AE = AH — q, where T is replaced by q in keeping 
with standard practice in percolation work and q is the 
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FIG. 11. The reduced pressure-probability phase diagram: 
the points show calculations performed at the probabilities 
below the critical point and the lines show fits to the Clau- 
sius-Clapeyron equation. 



average bond breaking probability considered; 0.56±0.02, 
0.54 ± 0.02 and 0.53 ± 0.02 for L = 9, 6 and 4 respec- 
tively. The AE/A values listed in Table VI were found 
by dividing AH — q by the most probable cluster size 



(1.15 ±0.05, 1.25 ±0.05 and 1.25 ±0.05 for L = 9, 6 and 
4 respectively), this puts AE on a "per site" basis. 

The values of AE/A shown in Table VI are nearly 
identical to the values of the surface energy coefficient 
Co , which is not surprising since for percolation on a sim- 
ple cubic lattice Co arises from the bonds broken to form 
the surface. Furthermore, the "energy of vaporization" is 
approximately equal to the number of bonds per lattice 



site (also shown in Table VI), a strong indication that the 
AE/A calculated here is the "bulk binding energy" of the 
lattice in question. The value of AE/A decreases with 
the size of the lattice because the percolation calculations 
were performed for open boundary conditions. 

The compressibility factor at the critical point was de- 
termined via Eq. (j|J) , the results are shown in Table jv|. 
From , q c and p c (determined below) the "pressure" 
at the critical point can be found. The resulting values 
of p c are shown in Table IV, but the interpretation of 
these values is an open question. 

Following the thermodynamic treatment of the perco- 
lation results, the reduced probability versus "reduced 
density" phase diagram is produced via Eq. (|3|). This 
leads to the points shown in Fig. [l2|. These points are 
then fit to Guggenheim's empirical formula, Eq. (32), 
with (3 — 0.43 ± 0.01 (in good agreement with text book 
values 0.41 |5^]) from Eq. (|3^). These results are shown 
for each lattice by the solid lines in Fig. H2. 




FIG. 12. The points shown on the low density branch are 
from the calculations performed at the bond breaking prob- 
abilities below the critical point. The lines are a fit to and 
reflection of Guggenheim's equation. The points shown on 
the high density branch show the size of the largest fragment 
at a given value of q normalized to the size of the largest 
fragment at q c . 



While it is not clear what density this plot describes, 
some insight can be gained by noting that the "liquid" 
branch reaches p/p c ~ 4 to ~ 4.5 at q = 0. Assuming 
that at q — p — 1 , since no bonds are broken, then p c ~ 
0.22 to ~ 0.25, which is approximately the percentage of 
bonds broken at the critical point. Thus it seems that 
the density in Fig. [l^ is related to the number of broken 
bonds. It is also noted that for q — the vapor branch 
of the coexistence curve shows p/p c > 0, this serves as 
an illustration of the magnitude of the error associated 
with this procedure. 

It is also possible to directly explore the behavior of the 
reduced density of the "liquid," at least in the larger sys- 
tem. This is done by normalizing the size of the largest 
cluster at a given value of q to the size of the largest 
cluster at the critical point A max (q) / A max (q c ) . Figure |l2| 
shows that for the L = 9 lattice, the measured normal- 
ized density of the liquid tracks along the coexistence 
curve predicted by Guggenheim's empirical formula and 
the reduced density of the vapor. For q/q c < 0.75 the 
effects of the finite size of the lattice are observed and 
the measured reduced density of the liquid deviates from 
the coexistence curve. The effects of finite size are more 
evident in the smaller lattices where there is little or no 
agreement between the measured reduced density of the 
liquid and the coexistence curves. Effects of finite size 
on the largest cluster such as these have been observed 
previously 59 1. 
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FIG. 13. The points are calculations performed at the bond 
breaking probabilities below the critical point and the lines 
are the results of the fits from the previous sections. 



For the sake of completeness, the "reduced pressure" 
versus "reduced density" projection of the phase diagram 
is shown in Fig. [13. 
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